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Abstract 

Einstein's field equations with variable gravitational and cosmological 
"constant" are considered in presence of perfect fiuid for Bianchi type-I 
space-time. Consequences of the four cases of the phenomenological decay 
of A have been discussed which are consistent with observations. The 
physical significance of the cosmological models have also been discussed. 
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1 Introduction 

There are significant observational evidence that the expansion of the Universe 
is undergoing a late time acceleration (Perlmutter et al. 1997, 1998, 1999; Riess 
et al. 1998, 2004; Efstathiou et al. 2002; Spergel et al. 2003; Allen et al. 2004; 
Sahni and Starobinsky 2000; Peebles and Ratra 2003; Padmanabhan 2003; Lima 
2004). This, in other words, amounts to saying that in the context of Einstein's 
general theory of relativity some sort of dark energy, constant or that varies only 
slowly with time and space dominates the current composition of cosmos. The 
origin and nature of such an accelerating field poses a completely open question. 
The main conclusion of these observations is that the expansion of the universe 
is accelerating. 

Among many possible alternatives, the simplest and most theoretically ap- 
pealing possibility for dark energy is the energy density stored on the vacuum 
state of all existing fields in the universe, i.e., py = where A is the cos- 

mological constant. However, a constant A cannot explain the huge difference 
between the cosmological constant inferred from observation and the vacuum 



1 



energy density resulting from quantum field theories. In an attempt to solve this 
problem, variable A was introduced such that A was large in the early universe 
and then decayed with evolution (Dolgov 1983). Cosmological scenarios with a 
time- varying A were proposed by several researchers. A number of models with 
different decay laws for the variation of cosmological term were investigated 
during last two decades (Chen and Hu 1991; Pavon 1991; Carvalho, Lima and 
Waga 1992; Lima and Maia 1994; Lima and Trodden 1996; Arbab and Abdel- 
Rahaman 1994; Vishwakarma 2001, Cunha and Santos 2004; Carneiro and Lima 
2005). 

On the other hand, numerous modifications of general relativity to allow for 
a variable G based on different arguments have been proposed (Wesson 1980). 
Variation of G has many interesting consequences in astrophysics. Canuto and 
Narlikar (1980) have shown that G-varying cosmology is consistent with what- 
soever cosmological observations available at present. A modification linking 
the variation of G with that of variable A-term has been considered within 
the framework of general relativity by a number of workers (Kallingas et al. 
1992; Abdel-Rahaman 1990; Herman 1991; Beesham 1986). This modification 
is appealing as it leaves the form of Einstein's equations formally unchanged by 
allowing a variation of G to be accompanied by a change in A. Cosmological 
models with time-dependent G and A in the solutions A A ~ were 

first obtained by Bertolami (1986). The cosmological models with variable G 
and A have been recently studied by several authors (Arbab 2003; Sistero 1991; 
Sattar and Vishwakarma 1997; Pradhan et al., 2001, 2002, 2005, 2007; Singh et 
al., 2006, 2007). 

Another important quantity which is supposed to be damped out in the 

course of cosmic evolution is the anisotropy of the cosmic expansion. The- 
oretical arguments and recent experimental data support the existence of an 
anisotropic phase that approaches an isotropic one. Therefore, it makes sense 
to consider the models of the universe with anisotropic background in the pres- 
ence of dark energy. 

The simplest of anisotropic models are Bianchi type-I homogeneous models 
whose spatial sections are flat but the expansion or contraction rate are direction 
dependent. For studying the possible effects of anisotropy in the early universe 
on present day observations many researchers (Huang 1990; Chimento et al. 
1997; Lima 1996; Lima and Carvalho 1994; Pradhan et al. 2004, 2006; Saha 
2005, 2006) have investigated Bianchi type-I models from different point of view. 

In the present article, we present a new class of solutions to Einstein's field 
equations with variable G and A in Bianchi type-I space-time in the presence of 
a perfect fluid. Consequences of the following four cases of the phenomenological 
decay of A have been discussed: 

Case 1 : A ~ iJ^ 
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Case 2 : A - iJ, 



Case 5 : A ^ p, 
Case 4 : A - R-^, 

where H, p, R are respectively the Hubble parameter, energy density and av- 
erage scale factor of the Bianchi type-I metric. The dynamical laws proposed 
for the decay of A have been widely studied by Chen and Wu (1990), Carvalho 
et al. (1992), Schutzhold (2002), Vishwakarma (2000), Arbab (1997, 1998) to 
name only a few. 



2 The Metric, Field Equations and Solutions 



We consider the space-time admitting Bianchi type-I group of motion in the 
form 

ds^ = -dt^ + A'^{t)dx^ + B^it)dy^ + C^{t)dz^. (1) 

We assume that the cosmic matter is represented by the energy-momentum 
tensor of a perfect fluid 



(2) 



where p, p are the energy density, thermodynamical pressure and Vi is the four- 
velocity vector of the fluid satisfying the relation 



v,v 



The Einstein's field equations with time-dependent G and A are 

1 



Ri 



-Rgi, 



-87rG{t)T,j + A{t)g,. 



(3) 



(4) 



For the metric ([T]) and energy-momentum tensor ^ in comoving system of 
coordinates, the field equation ^ yields 



(5) 

(6) 
(7) 
(8) 

(9) 
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8nGp - A 
SttGp + A = 

In view of vanishing divergence of Einstein tensor, we get 
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The usual energy conservation equation Tf., — yields 



, A B C 



(10) 



Equation © together with (fTU)) puts G and A in some sort of coupled field given 

by 

87rpG + A = 0. (11) 

Here and elsewhere a dot stands for ordinary time-derivative of the concerned 
quantity. From equation pT|) one concludes that when A is constant or A = 0, 
G turns out to be constant. 

Let R be the average scale factor of Bianchi type-I universe i.e. 



= = ABC. 

From equations ([5]), ([6]) and ([7]), we obtain 

A B _ ki 
A^B^R3' 

and 

B C _ k2 

b'c^r^' 

On integration equations (fT5|) and give 



A = mi R exp 

B — 7712 R exp 
C = 7773 exp 



(12) 

(13) 
(14) 
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dt 
R^ 

dt 
R^ 
d 

R3 



(15) 



where fci, ^2, 777i, 7772, 7773 are arbitrary constants of integration satisfying 

7)1177727773 — 1. 

Similar expressions as (fTS]) have also been established by Saha (2005). 

Hubble parameter _ff , volume expansion 9, shear a and deceleration param- 
eter q are given by 

e^3H = 3-^ 



V3R3' 

q = 



R' 

k > 0, (constant) 
H 



-1 - 
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Equations (HI)-® and pUj) can be written in terms of H, a and q as 

87rGp = i/2(2g-l)-(T2+A, (16) 
SttGp = Si/^ - 0-2 _ A, (17) 

P + 3(p + p)| = 0. (18) 

It is to note that energy density of the universe is a positive quantity. It is 
beUeved that at the early stages of the evolution when the average scale factor 
R was close to zero, the energy density of the universe was infinitely large. On 
the other hand, with the expansion of the universe i.e. with increase of i?, the 
energy density decreases and an infinitely large R corresponds to a p close to 
zero. In that case from ([T7|. we obtain 

3ij2 - A ^ 0. (19) 

From equation one concludes that: 

(i) A is essentially non-negative, 

(ii) in absence of a A-term beginning from some value of i?, the evolution of 
the universe becomes stand-still i.e. i? becomes constant since H becomes zero, 

(iii) in case of a positive A, the process of evolution of the universe never 
comes to a halt. Moreover, it is believed that the presence of dark energy (given 
by positive A) results in the accelerated expansion of the universe. As far as 
negative A is concerned, its presence imposes some restriction on p i.e p can 
never be small enough to be ignored. It means, in that case there exists some 
upper limit for i? as well. It is worth mention here that Saha (2006) has also 
given such conclusion in his paper but his approach was quite different. 
From equation (|17p . we obtain 

cr2 _ 1 SttGp a 
^ ^ 3 ^2 ffi- 

Therefore, < < ^ and < < i for A > 0. 

Thus, the presence of a positive A puts restriction on the upper limit of 
anisotropy whereas a negative A contributes to the anisotropy. 

From equation (|16p . we obtain 

Thus for negative A, the universe will always be in decelerating phase whereas 
a positive A will slow down the rate of decrease. Also (74 — —3aH implying 
that a decreases in an evolving universe and for infinitely large value of R, a 
becomes negligible. 
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Equations ([S]) - ([S]) and (fTT|) together with one of the decay laws for A 
given by cases (1) - (4) supply six equations in seven unknown functions of 
time A, B, C, p, p, A and G. To have deterministic solutions, we require one 
more condition. For this purpose, we assume that the volume expansion 9 is 
proportional to eigen values of shear tensor aij. It is believed that evolution 
of one parameter should also be responsible for the evolution of the others 
(Vishwakarma, 2005). Following Roy and Singh (1985), we take the volume 
expansion 9 having a constant ratio to the anisotropy in the direction of unit 
space-like vector A* i.e. ^ ^^^j is constant. In general, the above condition gives 
rise to 



(20) 



where m and n are constants. Using condition ((20|) in equations (fT3|) and (fT4|) . 
we obtain 



fcl-(m-l)fc2 
C = bl{k^t + fc4)('» + " + 2)'=l-('"-2"-l)'=2 



fci TO — 2n — 1 
J or — 7^ . 

k2 ni + n + 2 ' 



exp 
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= h exp 
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for 



ki TO — 2n — 1 
J or — f= . 

K2 TO + n + 2 ' 

;i TO — 2n — 1 



k2 



m + n + 2 



(21) 



(22) 



provided to + n 7^ 1. In the above ^3, k^, fcs and 61, 62, &3 are constants 
of integration. For these solutions, metric ^ takes the following forms after 
suitable transformations: 



2{m + n)ki+2nk2 2ki+2nk2 

TO — 2n — 1 



2fci-2(m-l)fc2 k^ 

_l_ 7-(™+„+2)f.i-(™-2„-i)fc2 dZ"^ for — 

k2 
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exp 



exp 



2^2 (to — n)T 



m + n + 2 
2k2(m+l)T' 



dX^ + exp 



TO + n + 2 



dZ^ for 
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(23) 



dY' 



m + n + 2 



(24) 



3 Discussion 



We now describe the models resulting from different dynamical laws for the de- 
cay of A. 
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For the model average scale factor R is given by 



Volume expansion 9, Hubble parameter H and shear a for the model are: 

1 

e = 3H=-, a2 = -^. 

T' 3T2 

Thus we see that j — Therefore, the model does not approach isotropy. If 
k is small, the models are quasi-isotropic i.e. ^ — 0. 



3.1 Case 1 

We consider 



A = 3PH^ 



where /3 is a constant of the order of unity. Here (3 represents the ratio between 
vacuum and critical densities. From equations ([5]), ([5]) and ([TT]) . we obtain 



87rp=^ ^^^^ (25) 

3fco 

Snp=^l^^±^T-^, (26) 
3fco 

A=^, (27) 

213 



G = koTi^-"''-!^) , fco > (constant). (28) 

We observe that the model has singularity at T = 0. It starts with a big bang 
from its singular state and continues to expand till T — oo. At T = 0, p, p, 
A, 9 and a are all infinite whereas G = for /3 > and G = oo for /? < 0. 
For infinitely large T, p, p, A, 9 and a are all zero but G = oo for (3 > and 
G = for P < 0. We also observe that in the absence of cosmological term 
A(/3 — 0), p — p i.e. matter content turns out to be a stuff fluid. For (3 > 0, 
p > p and p < p when /S < 0. When P — — 1, p — 0. The density parameter 
r2=^ = l — fc^ — /9 implying that Pc > P and pc < p for (3 > — fc^ and /3 < —k-^ 
respectively whereas pc = p when 13 = —k"^. The ratio between vacuum and 
matter densities is given by 

Pv_ _ 

p ~ l-k^~l3' 
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3.2 Case 2 : 

We now consider 



A = aH, 



where a is a positive constant of order of where m « IbQMeV is the energy 
scale of chiral phase transition of QCD (Borges and Carneiro, 2005). For this 
case, equations (O, (H]) and (fTTjl yield 

(1 - ef - a^T^ 

A^|., (31) 



The model has singularity at T = 0. The model starts from a big bang with p, 
p, A, 9, a all infinite and G finite. Thereafter p, p, A, 6 and a decrease and G 
increases. When T — we obtain p = 0, p = 0, A ^ilk'^) — ^i^-^-k'^) 

2 2 

and G is infinite. As T — + oo, p ^ „ ° , , p ^ — „ ° , , and 0, cr, G, A tend to 
zero. The density parameter fl = \ — k"^ — aT and the ratio between vacuum 
and critical densities is given by 



3.3 Case 3 : 

We now consider 



SnaGp 

A = 



3 ' 

where a is a constant. In this case from equations (O, ([8]) and (fTTj) . we obtain 

fl — k'^) 2(a + 3) 

8,p.(i_£l<|!i±lT-^, (34) 

ko(a + 6) 

^-3(a + 3)T2' ^"^^^ 

G^koT^. (36) 
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This model also starts from a big bang at T = with p, p, A, 9, a all infinite 
and G — (for a > 0) and it evolves to p ^ 0, p ^ 0, ^ 0, tr — *■ 0, A — » 
and G — > oo as T — > oo. The density parameter fl for this model is given by 



a + 3 ' 

and the ratio between vacuum and critical densities is obtained as 

py a{l - P) 



3.4 Case 4 : 

Finally we consider the case 



Pc a + 3 



A- ^ 



where 7 is a parameter to be determined from the observations. In this case 
from equations ([1]) - (I5|) and pT|) . we obtain 

s,,= Ki^i!)ri^. ,3.) 

8.p=Ki^i2Zli±M, (38) 

oka 

A = 7r-i, (39) 

G = ko[l-k^-3-fTi]^^. (40) 

Here we observe that this model also has singularity at T = 0. It starts from a 
big bang singularity with p, p, 9, A, a all infinite but G finite. For A > i.e. 

7 > 0, p becomes zero at T = (3:^) * whereas for A < i.e. 7 < 0, p = 

at T = (■^a^) . As T ^ 00, 9, a, A and G become zero but p and p become 
finite. The density parameter O for this model is given by 

Q. = l-k^ - 37T3. 

The ratio between vacuum and critical densities is given by 

^ = 37ri 

Pc 

The model (l24l) is not of much interest since it reduces to a static solution. 
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4 Conclusion 



In this paper, we have presented a class of solutions to Einstein's field equations 
with variable G and A in Bianchi type-I space-time in the presence of a perfect 
fluid. In some cases, it is observed that G is an increasing function of time. 
When the universe is required to have expanded from a finite minimum volume, 
the critical density assumption and conservation of energy-momentum tensor 
dictate that G increases in a perpetually expanding universe. The possibility of 
an increasing G has been suggested by several authors. 

We would like to mention here that Beesham (1994), Lima and Carvalho 
(1994), Kallingas et al. (1995) and Lima (1996) have also derived the Bianchi 
type I cosmological models with variable G and A assuming a particular form 
for G. These models have some similarities with our model ((23)) in the cases 
(1) and (3) only. But our derived results differ from these models in the sense 
that the both of these are constrained by the equation of state whereas we have 
neither assumed equation of state nor particular form of G. 

The behaviour of the universe in our models will be determined by the cos- 
mological term A ; this term has the same effect as a uniform mass density 
Peff = —K/AitG, which is constant in space and time. A positive value of A 
corresponds to a negative effective mass density (repulsion). Hence, we expect 
that in the universe with a positive value of A, the expansion will tend to ac- 
celerate; whereas in the universe with negative value of A, the expansion will 
slow down, stop and reverse. Recent cosmological observations (Garnavich et 
al. 1998; Perlmutter et al. 1997, 1998, 1999; Riess et al. 1998, 2004; Schmidt et 
al. 1998) suggest the existence of a positive cosmological constant A with the 
magnitude K{Gh/c') k, 10"^^'^. These observations on magnitude and red-shift 
of type la supernova suggest that our universe may be an accelerating one with 
induced cosmological density through the cosmological A-term. Thus, our mod- 
els are consistent with the results of recent observations. 
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